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fi]  The  coronal  index  (Cl)  of  solar  activity  is  the  irradiance  of  the  Sun  as  a  star  in  the 
coronal  green  line  (Fe  XIV,  530.3  nm  or  5303  A).  It  is  derived  from  ground-based 
observations  of  the  green  corona  made  by  the  network  of  coronal  stations  (currently 
Kislovodsk,  Lomnicky  Stit,  Norikura,  and  Sacramento  Peak).  The  Cl  was  introduced 
by  Rybansky  (1975)  to  facilitate  comparison  of  ground-based  green  line  measurements 
with  satellite-based  extreme  ultraviolet  and  soft  X-ray  observations.  The  Cl  since  1 965  is 
based  on  the  Lomnicky  Stit  photometric  scale;  the  Cl  was  extended  to  earlier  years  by 
Rybansky  et  al.  (1994)  based  on  cross-calibrations  of  Lomnicky  Stit  data  with 
measurements  made  at  Pic  du  Midi  and  Arosa.  The  resultant  1939-1992  Cl  had  the 
interesting  property  that  its  value  at  the  peak  of  the  11 -year  cycle  increased  more  or  less 
monotonically  from  cycle  1 8  through  cycle  22  even  though  the  peak  sunspot  number  of 
cycle  20  exhibited  a  significant  local  minimum  between  that  of  cycles  19  and  21. 

Rusin  and  Rybansky  (2002)  recently  showed  that  the  green  line  intensity  and  photospheric 
magnetic  field  strength  were  highly  correlated  from  1976  to  1999.  Since  the  photospheric 
magnetic  field  strength  is  highly  correlated  with  sunspot  number,  the  lack  of  close 
correspondence  between  the  sunspot  number  and  the  Cl  from  1939  to  2002  is  puzzling. 

Here  we  show  that  the  Cl  and  sunspot  number  are  highly  correlated  only  after  1 965, 
calling  the  previously-computed  coronal  index  for  earlier  years  (1939-1965)  into 
question.  We  can  use  the  correlation  between  the  Cl  and  sunspot  number  (also  the 
2800  MHz  radio  flux  and  the  cosmic  ray  intensity)  to  recompute  daily  values  of  the  Cl  for 
years  before  1966.  In  fact,  this  method  can  be  used  to  obtain  Cl  values  as  far  back  as  we 
have  reliable  sunspot  observations  (~1850).  The  net  result  of  this  exercise  is  a  Cl  that 
closely  tracks  the  sunspot  number  at  all  times.  We  can  use  the  sunspot-CI  relationship 
(for  1966-2002)  to  identify  which  coronal  stations  can  be  used  as  a  basis  for  the 
homogeneous  coronal  data  set  (HDS)  before  1966.  Thus  we  adopt  the  photometric  scale  of 
the  following  observatories  for  the  indicated  times:  Norikura  (1951-1954;  the  Norikura 
photometric  scale  was  also  used  from  1939  to  1954);  Pic  du  Midi  (1955-1959); 

Kislovodsk  (1960-1965).  Finally,  we  revised  the  post-1965  HDS  and  made  several  small 
corrections  and  now  include  data  from  Kislovodsk,  Norikura,  and  Sacramento  Peak  to  fill 
gaps  at  Lomnicky  Stit. 

Citation:  Rybansky,  M„  V.  Rusin,  M.  Minarovjech,  L.  Klocok,  and  E.  W.  Oliver  (2005),  Reexamination  of  the  coronal  index  of  solar 
activity,  J.  Geophys.  Res.,  110,  A08106,  doi: 1 0. 1 029/2005 JA0 11146. 


1.  Introduction 
1.1.  Background 

[2]  Young  and  Harkness  independently  made  the  first 
observation  of  the  coronal  emission  line  at  5303  A  during 
the  total  solar  eclipse  of  7  August  1869.  The  “green  line,” 
as  it  came  to  be  called,  was  originally  attributed  to  the 
hypothetical  element  “coronium.”  The  puzzle  of  the  5303 
A  line  and  coronium  was  solved  by  Grotrian  [1939]  and 
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Edlen  [1942]  who  concluded  that  the  emission  line  origi¬ 
nated  in  highly  ionized  Fe  (Fe  XIV). 

[3]  A  new  era  in  coronal  line  observations  began  in  1939 
when  Waldmeier  initiated  green  line  measurements  at  Arosa 
(Switzerland)  with  a  Lyot  [1932]  coronagraph  [Waldmeier, 
1957],  The  coronagraph  patrol  started  by  Waldmeier  was 
taken  up  over  the  years  at  other  coronal  stations  resulting  in 
an  inhomogeneous  set  of  records  that  extends  to  the  present 
time.  These  records  provide  an  indicator  of  solar  activity  in 
the  corona,  in  much  the  same  way  that  the  sunspot  record 
characterizes  photospheric  activity.  At  present,  four  coronal 
stations  (Lomnicky  Stit  (Slovakia),  Norikura  (Japan),  Kislo- 
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Figure  1.  Schematic  indicating  a  unit  volume  in  front  of 
the  solar  disk.  Integrating  the  volume  emissivity  over  the 
visible  disk  and  above  the  limb  of  the  Sun  yields  the  coronal 
index  (Cl). 


vodsk  (Russia),  and  Sacramento  Peak  (USA))  observe  the 
green  corona  regularly. 

[4]  The  green  coronal  line  is  the  most  prominent  of 
28  coronal  emission  lines  in  the  visible  part  of  the 
electromagnetic  spectrum.  It  is  observed  during  the  com¬ 
plete  11 -year  solar  cycle  activity  and  around  the  entire 
solar  limb.  Measurements  of  the  green  coronal  intensity 
represent  one  of  the  longest-running  direct  indices  of  solar 
activity,  being  surpassed  only  by  the  sunspot  number 
(available  since  1610),  sunspot  area  (1874),  and  calcium 
plage  index  (1905)  [Lean,  2000,  and  references  therein]. 
The  green  line  measurements  can  be  compared  with  these 
and  other  direct  solar  indices  (e.g.,  irradiance  observations 
from  space  [Frohlich  and  Lean,  1998;  Frohlich,  2003,  and 
references  therein])  as  well  as  with  indirect  indices  (e.g., 
cosmogenic  radionucleides  [Beer,  2000])  to  gain  insight  on 
solar  variability.  The  green  line  coronal  measurements  also 
have  potential  applications  for  understanding  the  impact  of 
solar  variability  on  climate  change  [e.g.,  Krivova  and 
Solanki,  2003,  and  references  therein],  coronal  hole  deter¬ 
minations  from  1939  and  for  improving  space  weather 
forecasting  [e.g.,  Mavromichalaki  et  al.,  2003]. 

1.2.  Coronal  Index  (Cl)  of  Solar  Activity 

[5]  To  facilitate  comparisons  of  ground-based  green  line 
measurements  with  observations  in  extreme  ultraviolet  and 
soft  X-ray  wavelengths  made  from  space,  Rybansky  [1975] 
proposed  a  new  index  of  solar  activity  based  on  Lomnicky 
Stit  530.3  nm  measurements  that  he  called  the  “coronal 
index  of  solar  activity”  (Cl).  The  Cl  expresses  the  irradi¬ 
ance  of  the  green  corona  for  the  Sun  as  a  star  and  can  be 
easily  compared  with  similar  full-disk  solar  indices.  Cl 
values  are  expressed  in  power  units  [W  sr~‘]  and  can  be 
easily  transferred  to  similar  units  obtained  from  satellite 
measurements  as  shown  below: 

10I6W  sr1  =4.5  x  ltr7W  m_2(at  I  AU) 

e  1.2  x  108photons  cm_2s_l.  (1) 


Monthly  averages  of  Cl  for  1939-2002  fluctuate  in  ranges 
from  2  x  1016  W  sr-1  (around  cycle  minimum)  to  20  x 
1016  W  sr-1  (around  cycle  maximum).  Daily  Cl  maximum 
values  have  never  exceeded  30  x  1016  W  sr-1.  The  green 
line  intensity  in  calibrated  observations  is  expressed  in 
absolute  coronal  units  (ACU).  One  ACU  represents  the 
intensity  of  the  continuous  spectrum  of  the  center  of  the  solar 
disk  in  the  width  of  1  A,  at  the  same  wavelength  as  the 
observed  coronal  spectral  line  (1  ACU  =  3.89  W  m~2  sr-1  at 
530.3  nm). 

[6]  According  to  an  agreement  between  coronal  observers 
in  1 947,  coronal  intensities  are  derived  at  least  once  per  day 
at  heights  of 40-60  arcsec  above  the  solar  surface  with  a  lag 
of  5  degrees  in  positional  angle.  The  initial  observation  is 
made  above  the  solar  north  pole,  and  71  subsequent  data 
points  are  made  at  the  same  height  proceeding  in  a  counter¬ 
clockwise  direction  (north,  east,  south,  west)  [d'Azambuja, 
1947].  The  Cl  for  a  given  day  is  based  on  limb  observations 
of  the  green  line  for  that  day  and  intervals  of  6  days  on  either 
side  of  that  day.  The  distribution  of  intensity  above  the  solar 
surface  is  obtained  by  using  the  observations  at  the  E-limb 
(W-limb)  from  the  6  days  preceding  (following)  the  given 
day  to  specify  the  coronal  intensities  above  the  eastern 
(western)  half  of  the  disk.  By  integrating  over  the  solar  disk 
(Figure  1),  we  obtain  the  irradiance  in  front  of  the  visible  part 
of  the  Sun  (EH)  to  which  we  add  the  irradiance  above  the 
solar  limb  (0.5  EL)  [see  Rybansky,  1 975,  Figure  2]  to  obtain 
the  coronal  index 

Cl  =  EH  +  0.5El  (2) 

which  is  the  total  irradiance  of  the  green  corona  into  one 
sterradian  (sr)  toward  Earth  [for  details  of  the  computation, 
see  Rybansky  [1975]. 

[7]  The  limitations  and  underlying  assumptions  of  the 
calculation  are  as  follows: 

[8]  1.  It  is  assumed  that  the  solar  corona  (green  line 
irradiance)  is  stable  for  the  6-day  period  required  for  regions 
to  pass  from  the  limb  to  central  meridian  (the  E-limb  case), 
or  from  central  meridian  to  the  limb  (W-limb  case). 

[9]  2.  Measurements  are  made  with  a  resolution  of 
5  degrees,  so  values  of  the  irradiance  at  intermediate 
points  must  be  interpolated. 

[10]  3.  The  existence  of  an  average  gradient  of  the  green 
coronal  intensity  is  assumed.  Thus  measurements  made  at 
the  height  of  40  arcsec  (as  in  the  case  of  Lomnicky  Stit) 
allow  us  to  derive  the  green  line  intensity  at  any  height 
above  the  solar  surface.  The  gradient  used  for  heights 
>40  arcsec, 

I(h)/I(40  arcsec)  =  0.473e~°'l28h  +  1 ,409e“"m  (3) 

was  taken  from  Rusin  [1973]. 

[11]  The  series  of  daily  coronal  intensities  dating  from 
1939  on  which  the  computation  of  Cl  is  based  is  referred 
to  as  the  homogeneous  coronal  data  set  (HDS).  For  each 
day,  the  HDS  consists  of  72  green  line  intensities  (cover¬ 
ing  all  azimuths)  which  have  been  transformed  to  the 
Lomnicky  Stit  photometric  scale  using  a  correlation  tech¬ 
nique  for  common  observational  days.  Missing  observa¬ 
tional  days  (short  gaps)  are  filled  in  by  interpolation  to 
give  a  more  continuous  data  set.  For  years  before  1966 
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(when  Lomnicky  Stit  observations  commenced),  cross¬ 
calibrations  were  used  in  turn  between  Lomnicky  Stit 
and  Pic  du  Midi  and  Pic  du  Midi  and  Arosa  to  extend 
the  HDS  back  to  1939.  See  Rybansky  et  al.  [1994,  2001] 
for  details  of  this  process. 

[12]  The  resulting  Cl  for  the  period  1939-2002 
exhibited  a  monotonic  increase  of  Cl  values  for  the  peak 
years  of  solar  cycles  18-22  (1944  to  1996),  while  the 
corresponding  peak  sunspot  numbers  had  a  local  minimum 
for  solar  cycle  20.  Recently,  Rusin  and  Rybansky  [2002] 
showed  that  the  green  line  intensity  and  photospheric 
magnetic  field  strength  were  highly  correlated  from  1976 
to  1999.  Since  the  photospheric  magnetic  field  strength  is 
highly  correlated  with  sunspot  number,  the  lack  of  close 
correspondence  between  the  sunspot  number  and  the  Cl 
from  1939  to  1965  calls  the  pre-1966  Cl  data  into 
question. 

1.3.  Outline  of  Analysis 

[n]  In  this  paper,  we  use  correlations  between  the  Cl  and 
three  other  indices,  sunspot  number,  2800  MHz  radio  flux, 
and  cosmic  ray  intensity,  for  the  years  1966-2002,  to 
recreate  the  Cl  for  earlier  years.  We  then  use  the  pre-1966 
Cl  values  to  determine  which  coronal  stations  to  use  as  a 
baseline  for  the  HDS  for  years  from  1939  to  1965.  Finally, 
we  use  a  more  comprehensive  data  base  to  compute  the 
HDS  after  1965  and  apply  some  additional  small  corrections 
to  the  1966-2002  HDS.  The  principal  goal,  and  end  result, 
of  this  work  is  a  revised  and  updated  HDS  and  Cl  series  for 
use  by  solar,  and  solar-terrestrial,  physicists. 

2.  Analysis 

2.1.  Coronal  Green  Line  Data 

[14]  Measurements  of  coronal  intensities  from  the  early 
coronagraph  stations  at  Arosa  [  Waldmeier ,  1957]  and  Pic  du 
Midi  [e.g.,  Leroy  and  Trellis,  1974]  differed  in  their  partic¬ 
ulars,  e.g.,  method  of  photometry  (visual  or  photographic, 
spectral  or  narrow  passband  filters),  calibration  (artificial 
lamp,  personal  estimation,  or  comparison  to  the  center  of 
the  solar  disk),  and  height  (above  the  limb)  of  observation. 
These  various  differences  or  inhomogeneities  between 
green  line  patrols  persisted  as  new  stations  came  on  line, 
even  when  efforts  were  made  to  standardize  observational 
methods  [d’Azambuja,  1947]. 

[15]  The  main  differences  between  stations  today  involve 
calibration  and  the  observational  method  used.  For  example, 
while  some  coronal  stations,  e.g.,  Norikura,  Pic  du  Midi, 
Climax,  Kislovodsk,  related  coronal  measurements  to  the 
intensity  observed  at  the  center  of  the  solar  disk,  other 
stations,  e.g.,  Arosa,  Wendelstein,  Kanzelhohe,  Climax 
(early  years  from  1942  to  1955),  Sacramento  Peak  (early 
years  from  1955  to  1966),  made  only  qualitative  visual 
observations  and  reported  intensities  on  an  arbitrary  scale 
from  0  to  50  with  only  occasional  calibration.  Intensities 
obtained  with  this  latter  method  are  called  “relative  inten¬ 
sities.”  To  create  the  HDS,  the  relative  intensities  must  be 
transformed  to  the  absolute  coronal  units  (ACU)  described 
in  section  1.  However,  this  is  difficult  to  do  in  practice  as 
already  noted  by  Bruzek  [1955],  Thus  in  constructing  the 
HDS,  comparatively  little  use  was  made  of  the  relative 
intensity  data. 


[16]  Even  the  absolute  coronal  intensities,  however, 
have  their  limitations  since  the  values  obtained  at  differ¬ 
ent  observatories  depend  on:  (1)  the  accuracy  of  the 
determination  of  the  optical  density  of  the  calibration 
filter  (the  optical  density  of  the  calibration  filter  at 
Lomnicky  Stit  is  ~4,  indicating  a  104  reduction  in  light 
from  the  center  of  the  Sun,  but  in  practice,  this  parameter 
is  difficult  to  determine  precisely  because  of  the  compli¬ 
cated  light  path  in  the  coronagraph),  and  (2)  the  reduction 
of  the  scattered  photospheric  light  (the  continuum  and  the 
blended  absorption  at  530.23  nm)  in  Earth’s  atmosphere 
and  the  instrument. 

[17]  The  accuracy  of  coronal  intensity  measurements 
depends  on  the  observational  height  above  the  solar  photo¬ 
sphere  because  of  the  gradient  of  the  green  corona  intensity 
(see  Rusin  [1973]  for  details).  A  change  of  the  observed 
height  above  the  solar  limb  by  1  arcmin  results  in  an 
intensity  decrease  by  approximately  a  factor  of  two.  (For 
a  small  (e.g.,  F  =  1.5  m)  coronagraph,  the  uncertainty  in 
pointing  accuracy  can  be  ~1  arcmin.)  Another  factor  that 
influences  the  coronal  intensity  is  the  device  used  to  isolate 
a  narrow  portion  of  the  spectrum:  (1)  a  spectrograph  or  (2)  a 
narrow  passband  filter.  It  is  generally  believed  (known)  that 
spectral  observations  calibrated  to  the  center  of  the  solar 
disk  using  a  spectrograph  have  higher  accuracy  than  those 
made  with  a  narrow  passband  filter.  The  above  reasons, 
coupled  with  real  changes  in  the  emission  corona  irradiance, 
which,  as  we  know,  may  occur  on  short  timescales  (relative 
to  the  difference  in  time  at  which  stations  spanning  the 
Earth  make  their  observations),  cause  coronal  intensities 
from  different  coronal  stations  to  vary.  An  indication  of  the 
extent  of  the  discrepancy  can  be  seen  in  Figure  2  where  we 
plot  highly-smoothed  (81 -day  running  averages)  data  for  all 
of  the  coronal  stations  considered  in  this  study,  for  the 
intervals  during  which  they  observed  between  1939  and 
2002.  The  large  station-to-station  differences  evident  can 
not  be  explained  in  terms  of  the  time  shift  between  indi¬ 
vidual  coronal  stations;  any  discrepancies  due  to  differences 
in  the  UT  time  of  observation  will  be  greatly  reduced  by  in 
the  long-term  averages. 

2.2.  A  Difference  in  Values  of  the  Coronal  Index 
Before  and  After  1966.0 

[is]  In  order  to  examine  the  stability  of  the  Lomnicky  Stit 
photometric  scale  after  1965,  and  to  check  if  earlier  pre- 
1966  Cl  data  published  by  Rybansky  et  al.  [1994]  are 
correct,  we  did  a  correlation  between  monthly  averages  of 
the  Cl  (Cl),  sunspot  number  (R),  the  2800  MHz  radio  flux 
(RF),  and  cosmic  ray  (CR)  intensity  for  all  available 
common  years  and  also  for  the  identical  time  periods 
1953-2002  as  shown  in  Table  1.  The  empirical  relations 
between  the  Cl  and  the  three  other  parameters  are  given  in 
Table  1,  along  with  the  corresponding  correlation  coeffi¬ 
cients.  In  all  cases,  the  correlation  coefficients  are  relatively 
high. 

[19]  Using  the  relations  in  Table  1  we  can  compute  a  new 
series  of  CIs,  plotted  as  dotted  lines  in  Figure  3  and  marked 
as  Cl  (R),  Cl  (RF),  and  Cl  (CR)  for  comparison  with  the 
current  index  (dashed  lines).  Inspection  of  Figure  3  reveals 
that  the  correlation  is  much  weaker  for  cycles  before  1966 
than  after.  This  behavior  is  another  proof  that  the  Lomnicky 
Stit  photometric  scale  in  the  period  1966-2002  is  stable, 
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Figure  2.  The  81 -day  running  averages  of  the  green 
coronal  line  intensities  for  different  coronal  stations, 
showing  differences  in  absolute  values  and  trends  attribu¬ 
table  to  various  factors  discussed  in  the  text. 

and,  on  the  opposite  hand,  it  suggests  a  problem  with  the 
pre-1966  Cl  values. 

3.  HDS  and  Cl  Reconstruction  Before  1966 

[20]  Given  the  results  in  Figure  3,  it  is  clear  that  the 
current  HDS,  on  which  Cl  is  based,  needs  to  be  modified  for 
years  before  1966.  The  current  pre-1966  HDS  uses  data 
from  Pic  du  Midi  and  Arosa  as  reference  stations,  with  their 
data  adjusted  to  the  Lomnicky  Stit  photometric  scale.  The 
basic  questions  in  this  reanalysis  are  as  follows:  (1)  How  to 
select  observational  days  from  the  data  sets  for  the  individ¬ 
ual  coronal  stations?  (2)  Which  subset  of  coronal  stations 


Table  1.  Relations  Between  Monthly  Averages  of  Cl  and  R,  RF, 
and  CR° 


Relation 

Period 

CC 

CI(R)  =  3.64  +  0.0548  x  R 

1939-  2002 

r  =  0.729 

=  3.31  +  0.0611  x  R 

1953-2002 

r  =  0.782 

Cl  (RF)  =  -0.219  +  0.00603  x  RF 

1947  -  2002 

r  =  0.754 

=  -0.766  +  0.00667  x  RF 

1953-2002 

r  =  0.796 

Cl  (CR)  =  59.6  —  0.0131  x  CR 

1953-2002 

r  =  0.813 

“Correlation  coefficients  (CC)  are  shown  in  the  last  column. 


(Pic  du  Midi,  Arosa,  Sacramento  Peak,  Climax,  Norikura, 
Wendelstein,  Kanzelhohe,  Kislovodsk)  should  be  used  as 
reference  stations  during  the  pre-1966  period?  (3)  How  to 
transform  the  photometric  reference  scale(s)  of  the  new 
reference  station(s)  to  the  Lomnicky  Stit  scale?  These  steps 
are  required  to  create  a  reconstructed  homogeneous  coronal 
data  set  over  the  entire  1939-2002  period,  and  to  recom¬ 
pute  the  pre-1966  Cl.  Our  procedure  is  given  below. 

3.1.  Data  Selection 

[21]  In  this  study,  we  considered  all  available  coronal 
green  line  data  (the  530.3  nm  emission  line  corresponds  to  a 
coronal  temperature  of  ~2  million  K)  from  the  various 
observatories  that  have  made  such  measurements  over  the 
years  as  published  in  the  Quarterly  Bulletin  of  Solar 
Activity.  However,  we  did  not  use  any  data  from  the  Ulan 
Bator  (Mongolia)  and  Alma  Ata  (Kazakhstan)  which  are 
generally  regarded  as  to  be  of  low  accuracy.  To  assist  in  the 
decision  about  which  data  (coronal  intensities)  to  consider 
from  the  remaining  coronal  stations,  we  created  an  auxiliary 
parameter  AI  for  each  observational  day  for  each  station.  AI 
is  the  daily  intensity  averaged  over  all  72  positional  angles. 
We  eliminated  observational  days  which  exhibited  large 
day-to-day  jumps  (~50%  or  more)  in  AI  values  or  for  which 
only  partial  data  (e.g.,  covering  only  half  of  the  positional 
angles)  were  available. 

[22]  After  removing  low  quality  data  days,  the  basic 
green  line  data  sets  having  intensities  in  absolute  coronal 


Figure  3.  Monthly  averages  of  Cl  from  1939  to  2002 
(heavy  lines  in  each  of  the  panels)  compared  with  proxy  CIs 
calculated  from  the  relationships  given  in  Table  1  (thin  lines 
in  each  of  the  panels). 
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Figure  4.  Correlation  coefficients  for  cross-correlations 
between  AI  and  R  for  a  range  of  lag  times  (9-days  running 
mean  is  shown  by  the  upper  curve,  3-days  running  mean  is 
shown  by  the  middle  curve;  and  1  day  running  mean  is 
shown  by  the  lowest  curve). 

units  that  were  used  to  reconstruct  the  1939-1965  HDS  and 
recomputed  the  Cl  since  1965  were  as  follows: 

[23]  1.  Pic  du  Midi  (PdM),  1943-1974;  we  used  only 
2992  selected  days  (from  3493  total  observational  days). 
PdM  changed  its  method  of  photometry  twice  during  this 
30-year  interval  [Leroy  and  Trellis ,  1974],  A  calibration 
problem  at  PdM  involving  an  east-west  asymmetry  (higher 
west  limb  intensities,  not  confirmed  by  other  stations)  was 
discussed  by  Pajdusakova  [1966]. 

[24]  2.  Norikura  (NO),  1951-1965;  only  1052  (of  1540) 
observational  days  are  used  in  this  period.  Later  data  (post- 
1965),  as  shown  by  Sakurai  et  al.  [1999],  were  not  as  stable 
in  absolute  units,  and  data  from  this  period  were  used  less 
frequently  (~42%  of  2864  observational  days). 

[25]  3.  Kislovodsk  (KI),  1957-2002;  4788  (of  5273) 
observational  days. 

[26]  4.  Climax  (CX),  1955-1965;  1427  (of  1930)  days. 

[27]  5.  Sacramento  Peak  (SP),  1958-1966;  913  (of  1218) 
days. 

[28]  Relative  (noncalibrated)  data  were  taken  from  fol¬ 
lowing  coronal  stations:  (1)  Arosa  (AR),  1939-1974;  2023 
(of  2340)  observational  days,  (2)  Climax  (CX),  1947— 
1954;  1371  (of  1852)  observational  days.  (3)  Wendelstein 
(WE),  1947-1978;  1200  (of  1467)  observational  days. 
(4)  Kanzelhohe  (KH),  1947-1963;  634  (of  812)  observa¬ 
tional  days.  The  relatively  low  number  of  observational 
days  at  each  station  for  the  intervals  considered  primarily 
reflects  the  rarity  of  sky  conditions  suitable  for  coronal 
observations. 

3.2.  Selection  of  Reference  Stations  for  the  Pre-1966 
Epoch 

[29]  In  Figure  4,  we  have  plotted  the  correlation  coeffi¬ 
cients  for  cross-correlations  between  sunspot  number  R  and 
AI  computed  from  the  old  HDS  for  1966-2002  for  all 
13,514  days  for  a  range  of  lead/lag  times  and  three  different 
data  averaging  intervals  (9  days,  3  days,  and  1  day).  The 
relation  between  these  two  parameters  is:  AI(R)  =  0.2R  + 


8  with  the  correlation  coefficient  0.737  when  the  daily  AI  is 
compared  with  R  for  7  days  later.  Using  this  relation,  we 
can  obtain  an  average  AI(R)  from  sunspot  number  alone. 

[30]  A  comparison  of  AI  values  for  Kislovodsk,  Norikura, 
and  Pic  du  Midi  along  with  AI(R)  for  1939-1966  is  shown 
in  Figure  5.  The  main  criterion  for  a  choice  of  reference 
station  for  the  pre-1966  epoch  is  a  good  correlation  with 
AI(R).  Thus  on  the  basis  of  the  data  in  Figure  5,  we  selected 
the  following  reference  stations  for  the  pre-1966  period: 
Kislovodsk:  1960-1965,  Pic  du  Midi:  1955-1959,  Nori¬ 
kura:  1951-1954. 

[31]  For  1939-1950,  data  from  Pic  du  Midi,  Climax, 
Arosa,  Wendelstein,  and  Kanzelhohe  were  converted  to  the 
Norikura  photometric  scale  using  correlations  from  com¬ 
mon  observational  days  between  Norikura  and  Pic  du  Midi 
in  1 95 1  —  1 954.  Rescaled  Pic  du  Midi  data  then  were  used  as 
reference  data  in  the  period  1943-1950.  Arosa  data  were 
transformed  to  the  Pic  du  Midi  rescaled  photometric  scale  in 
the  period  1939-1942. 

3.3.  Putting  the  Pre-1966  Data  on  the  Lomnicky  Sti't 
Photometric  Scale 

[32]  Different  coronal  intensities  from  individual  coronal 
stations  cannot  be  used  directly  to  reconstruct  the  HDS  and 
compute  the  Cl  [see  Rybansky,  1975]  but  first  must  be 
transformed  to  a  single  photometric  scale.  Based  on  the 
results  of  Figure  3,  we  use,  as  was  done  previously 
[Rybansky  et  al.,  1994],  the  Lomnicky  Stit  photometric 
scale  for  years  after  1965  to  recompute  Cl.  For  earlier 
years,  we  used  the  photometric  scales  of  the  reference 
stations  listed  in  section  3.2. 

[33]  A  further  comment  on  the  photometric  rescaling 
process  is  in  order.  As  a  refinement  on  our  process,  we 
took  one  additional  step  before  rescaling  the  data.  It  is  very 
easily  recognized  (see  Figure  6)  that  coronal  intensities 
between  individual  coronal  stations  differ  not  only  in  their 
intensities,  but  in  the  position  of  maximum  intensity  as  well 


Figure  5.  The  81-day  running  averages  (heavy  lines)  of 
AI  for  Kislovodsk  (1960-1965);  Pic  du  Midi  (1955-1959); 
Norikura  (1951-1954),  rescaled  Pic  du  Midi  (1944-1950), 
rescaled  Arosa  (1939-1943)  and  AI(R)  for  the  entire  period 
is  in  the  thin  line.  See  text  for  details. 


5  of  9 


A08106 


RYBANSKY  ET  AL.:  CORONAL  INDEX  OF  SOLAR  ACTIVITY 


A08106 


N 


S 


Figure  6.  Distribution  of  the  green  coronal  intensities  for 
common  observational  days  from  LS  (heavy  line)  and  KI 
(dotted  line  before  the  shift,  thin  line  after  the  10  degree 
shift). 

that  can  differ  for  E-  and  W-  limbs  [e.g.,  Gnevyshev,  1963; 
Sykora,  1971;  Rybansky,  1975].  Thus  we  correlated  the 
shapes  of  the  corona  for  common  observational  days  for 
one-year  periods  (because  the  angular  shifts  may  vary  over 
time)  to  determine  if  the  data  from  other  stations  needed  to 
be  rotated  to  conform  to  Lomnicky  Stft  (or  other  reference 
station  for  pre-1966  periods)  observations. 

3.4.  An  Example  of  the  Photometric  Rescaling  Process 

3.4.1.  Deciding  Which  Station  to  Use 

[34]  For  days  when  data  are  not  available  from  the 
reference  station,  it  is  necessary  to  decide  which  station  to 
use  to  fill  in  the  gap,  assuming  data  from  more  than  one 
stations  are  available.  The  choice  is  based  on  the  best 
correlation  between  reference  and  other  stations  (if  several 


Figure  7.  Comparison  of  green  line  coronal  intensities 
between  Lomnicky  Stft  (y  axis)  and  Kislovodsk  (x  axis)  for 
the  period  1966-1974.  Dots  indicate  individual  intensities 
for  each  positional  angle  for  common  observational  days; 
crosses  represent  averaged  intensities  (AI)  from  common 
observational  days.  The  thick  regression  line  is  based  on  a 
weighted  LSQ  and  the  dashed  lines  represent  99% 
confidence  levels. 


Figure  8.  Plot  of  81 -day  running  averages  of  AI  as 
derived  from  former  HDS0id  is  shown  by  the  thin  line 
(upper  panel),  revised  HDSrev  is  shown  by  the  thick  line 
(both  panels)  and  AI(R)  as  derived  from  sunspot  number  is 
shown  by  the  thin  line  (bottom  panel). 

stations  have  observations  for  the  some  day)  for  common 
observational  days  over  short  (~1  year)  intervals  including 
the  day  of  interest.  If,  however,  only  data  from  one  station 
are  available,  these  data  are  used. 

3.4.2.  Angular  Adjustment 

[35]  The  averages  from  common  observational  days  at 
each  positional  angle  for  the  reference  and  adopted  stations 
are  the  basis  for  the  shift  in  the  positional  angle.  Usually,  a 
polar  diagram  (Figure  6)  is  prepared,  and  the  data  for  the 
adopted  station  are  rotated  until  maxima  agree  with  that  of 
the  reference  station  (using  interpolation  as  needed). 

3.4.3.  Intensity  Scaling 

[36]  The  method  for  rescaling  the  intensity  from  the 
chosen  station  to  the  reference  photometric  station  is  illus¬ 
trated  in  Figure  7.  The  symbols  used  in  Figure  7  are 
interpreted  as  follows:  small  dots,  all  measured  intensities 
of  the  green  corona  for  common  observational  days  and  in 
all  positional  angles  in  1966-1974;  crosses,  averaged 
intensities  (AI)  from  common  observational  days  for  chosen 
coronal  station. 

[37]  To  prepare  a  conversion  relation,  indicated  by  the 
heavy  line  in  Figure  7,  we  used  a  simple  weighted  least 
squares  (LSQ)  straight  line  prepared  from  the  dots,  where 
the  weight  of  individual  dots  is  dependent  on  the  perpen¬ 
dicular  distance  from  a  LSQ  straight  line  (not  shown) 
derived  from  the  crosses. 

[38]  For  example,  the  conversion  relation  (derived  from 
6048  points  -  84  common  observational  days  for  1966— 

1 974)  between  LS  and  KI  is 

LS  =  0.62  KI  +  0.28  (4) 

There  were  557  LS  and  956  KI  observational  days  during 
this  9  year  interval  and  a  total  of  147  common  observational 
days.  However,  KI  observations  were  complete  (72  values) 
for  only  84  of  the  common  days.  This  process  was  done  for 
all  reference  coronal  stations  for  which  coronal  intensities 
were  transformed  from  other  stations. 

[39]  We  believe  that  above  method  can  be  used  to  accu¬ 
rately  convert  coronal  intensities  from  the  worldwide  net- 


6  of  9 


A08106 


RYBANSKY  ETAL.:  CORONAL  INDEX  OF  SOLAR  ACTIVITY 


A08106 


Figure  9.  The  total  number  of  interpolated  days  (thick 
line)  and  maximum  continuous  interpolated  interval  (thin 
line)  in  individual  years.  There  are  no  observational  data  in 
interpolated  periods. 

work  of  coronal  stations  to  the  Lomnicky  Stit  photometric 
scale,  and  that  the  resulting  homogeneous  coronal  data  set 
and  the  Cl  values  computed  from  it  are  now  correct.  A  plot 
of  81 -day  running  averages  of  AI  as  derived  from  the  daily 
old  HDS  means,  the  daily  revised  HDS  means,  and  AI(R)  are 
shown  in  Figure  8. 

3.5.  Filling  in  the  Data  Gaps 

3.5.1.  Additional  Data  for  the  Post-1965  Epoch 

[40]  To  make  the  HDS  as  complete  as  possible,  we 
have  added  additional  data  sets  since  the  initial  work  by 
Rybansky  et  al.  [1994],  The  new  data  sets  include  recent 


data  from  Kislovodsk  (1966-2002),  Norikura  (1966-2001) 
(T.  Sakurai,  private  communication,  2002);  and  Sacramento 
Peak  (1974-2002;  Solar  Geophysical  Data,  Prompt 
Reports).  These  data  are  folded  into  the  HDS  using  the 
photometric  rescaling  process  described  in  sections  3.3  and 
3.4.  The  additional  data  reduced  the  number  and  length  of 
intervals  with  no  data  which  were  previously  filled  by 
interpolation.  In  addition,  using  more  coronal  stations 
enabled  us  to  easily  recognize  and  eliminate  unreliable 
data. 

3.5.2.  Interpolation  Across  Data  Gaps 

[<u]  There  are  23,376  days  in  the  interval  from  1  January 
1939  and  31  December  2002.  The  total  number  of  days  with 
useful  green  coronal  data  used  to  create  the  original  HDS 
was  14,824  (63%  of  23,376)  versus  15,506  (66%)  in  the 
revised  HDS.  There  are  20,454  days  in  the  period  1947— 
2002,  and  14,898  (73%)  observational  days  were  used  in  the 
revised  HDS.  Coronal  intensities  for  missing  days  were 
filled  with  a  linear  interpolation  based  on  the  last  and  first 
day  of  observations  around  the  interpolated  period  (missing 
day  or  days). 

[42]  The  total  number  of  interpolated  days  (thick  line)  and 
the  maximum  interval  in  days  of  continuously  interpolated 
data  (thin  line)  in  individual  years  are  shown  in  Figure  9. 
After  1947  the  median  number  of  interpolated  days  per 
year  is  ~100  and  the  longest  continuous  interpolated 
interval  is  21  days.  Interpolated  intensities  for  missing 
days  are  indicated  by  using  a  “0”  for  the  coronal  station 
code  in  the  HDS. 

3.6.  Recomputed  Coronal  Index 

[43]  With  the  updated  HDS,  we  recomputed  Cl  for  the 
1939-2002  interval.  The  recomputed  Cl  (CIrev)  is  shown  in 


Figure  10.  Plot  of  Cl  monthly  means  (old  and  revised)  for  1939-2002  (top).  Difference  between 
CIoM  and  CIrev  (bottom). 
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Figure  11.  Correlation  between  Cl  and  the  sunspot 
number,  (left)  Correlation  between  the  CIoW  and  the  sunspot 
number  for  1939-1965  (dots  and  thin  line);  circles  and 
heavy  line  for  1966-2002.  (right)  The  same  for  CIrev  and 
the  sunspot  number. 

Figure  10  along  with  the  old  Cl  (C[0)a).  In  addition  to 
removing  the  pre-1966  discrepancies  with  the  sunspot 
number  (Figure  3),  we  see  that  incorporation  of  the  new 
data  also  results  in  improvement  in  the  post- 1965  period, 
particularly  for  the  1979-1983  interval.  CI0id  (left)  and 
CIrev  (right)  values  are  plotted  against  sunspot  number  for 
periods  before  and  after  1965  in  Figure  11. 


[44]  The  relation  between  the  revised  Cl  and  the  sunspot 
number  (for  81 -day  means)  is 

Cl  =  0.0837  R  +  2.04  (5) 

with  a  correlation  coefficient  of  0.914  in  the  period  1939- 
2002,  which  compares  with  a  correlation  coefficient  of 
0.729  for  a  correlation  of  the  sunspot  number  and  the  old 
Cl.  Using  the  relation  in  equation  (5),  the  Cl  can  be 
extended  to  1848  as  shown  in  Figure  12.  Note  that  because 
the  sunspot  number  returns  to  near  zero  values  at  each 
sunspot  minimum  that  the  reconstructed  Cl  has  a  nearly 
constant  floor  at  ~2.04.  The  relation  between  the  revised  Cl 
and  the  2800  MHz  radio  flux  (for  81 -day  means)  is 

Cl  =  0.00874  RF  -  3.06  (6) 

with  a  correlation  coefficient  0.924  (slightly  higher  than  for 
the  sunspot  number)  in  the  period  1947-2002. 

4.  Results  and  Discussion 

[45]  This  study  was  motivated  by  the  fact  that  the 
previously  computed  coronal  index  (Cl)  of  solar  activity 
[Rybansky  et  al.,  1994]  did  not  agree  well  with  the 
2800  MHz  radio  flux,  the  cosmic  rays  and  the  sunspot 
number  for  years  before  1966.  This  runs  counter  to  the 
recent  finding  by  Rusin  and  Rybansky  [2002]  that  the  green 
line  intensity  and  photospheric  magnetic  field  strength  were 
highly  correlated  from  1976  to  1999.  To  adjust  the  homo¬ 
geneous  coronal  data  set  on  which  Cl  is  based,  we  used 


Year 

Figure  12.  Monthly  averages  of  Cl  (1848-2002).  Heavy  line  shows  CIrev  (1939-2002)  and  thin  line 
shows  Cl  extrapolated  from  the  sunspot  number  (1848-1938). 
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correlations  between  coronal  intensities  and  the  sunspot 
number  to  identify  suitable  reference  stations  (yielding  better 
agreement  between  computed  Cl  and  sunspot  number)  for 
years  before  1966.  We  retract  our  previous  statement  [Rusin 
et  al,  2004],  based  on  the  former  Cl  data,  that  the  Cl  (at  solar 
maximum)  rose  steadily  between  1947  and  1996.  Neverthe¬ 
less,  like  the  sunspot  number,  Cl  can  vary  from  cycle  to  cycle 
by  ~30  percent  with  irregular  trends  for  longer  periods,  e.g., 
the  increase  in  the  period  1965-1996  (cycles  20,  21,  and 
22).  The  revised  and  extended  Cl  has  applications  for 
correlative  studies  such  as  those  of  Mavromichalaki  et  al. 
[2002,  2003], 

[46]  Our  decision  to  modify  the  pre-1966  HDS  (and  Cl  as 
a  result)  is  based  on  practical  experience  with  the  observa¬ 
tions  and  awareness  of  the  various  causes  of  discrepancies 
between  stations  and  the  recognition  that  the  modem  (post- 
1965)  data  at  Lomnicky  Stit  tracks  the  standard  measures  of 
photospheric  activity  (sunspot  number  and  10-cm  radio 
flux)  and  also  the  inverted  cosmic  ray  modulation  curve 
quite  well.  This  leads  us  to  believe  that  the  modern 
measurements  at  Lomnicky  Stit  more  accurately  reflect 
the  Sun’s  behavior,  and  we  have  normalized  the  earlier 
observations  accordingly.  Coronal  intensities  from  world¬ 
wide  stations  used  in  this  analysis  are  available  at  http:// 
www.ta3.sk/data/corona/index  and  http://www.noaa.edu/ 
corona,  as  are  the  reconstructed  HDS  (1939-2002)  and 
the  Cl  series  derived  from  it. 
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